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Abstract 

A phenomenological description of the high- Tc superconductors based on the 
Cooperon propagator is presented. This model ahows one to study the effects of 
local pairing correlations and long-range phase fluctuations on the same footing, 
both above and below Tc. Based on numerical calculations, it is shown that the two 
types of correlations contribute to the gap/pseudogap in the single-particle excita- 
tion spectra. The concourse of these two effects can induce low energy states, which 
should be observable in underdoped materials at very low temperature. 
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The nature of the pseudogap phase in the high-Tc cuprate superconductors (HTS) 
is still the subject of lively debate and its elucidation may be one of the keys to the 
understanding of these materials. At present we may perhaps, from the experimen- 
tal point of view, state the following: (1) The pseudogap — defined as a depletion of 
the single particle density of states (DOS) around the Fermi energy, which appears 
gradually at a temperature T* — goes smoothly into the superconducting gap at 
Tc, as seen in STM[1] and ARPES[2] measurements. (2) The STM measurements 
of vortices below Tc show that the DOS in the center of the vortex core ressembles 
strongly the DOS observed above Tc in the pseudogap regime[3,4]. (3) The fluctua- 
tions above Tc (and outside the 3D critical regime near Tc) seem to correspond to a 
phase fluctuation regime, quantitatively well described by the Kosterlitz-Thouless 
(KT) theory above Tkt[5,6]. (4) This phase fluctuation regime seems to disappear 
at a characteristic temperature well below T*[5,6]. One is therefore in need of a 
theoretical framework which (i) is able to show a smooth transition of the DOS 
across Tc, (ii) connects the phase fluctuation regime with the properties of the KT 
theory, and {in) relates the properties within a vortex core to the properties of the 
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pseudo-gapped phase. Such a theoretical framework has been developed 40 years 
ago by Kadanoff and Martin (KM) [7], who presented a theory of superconductivity 
entirely based on the pair correlation function (Cooperon propagator). In a recent 
paper [8], we applied this framework for the s-wave symmetry. The purpose of this 
communication is to extend these results to the d-wave symmetry case. 

In the KM formalism, the relation between the two-body Cooperon propagator 
L(r, r'; s', s; r) = (rT-{V'|(r,r)'i/'|(r',r)V'J(s,0)V'|(s',0)}) and the self-energy is 

S(r, s, r) = - / dr'ds' V{r, r')L{r, r'; s', s; t)V{s', s) Go(s', r', -r), (1) 



where Go is the free Green's function and V{r, r') is the effective interaction between 
the electrons. In a translationally invariant system with a short-ranged interaction, 
we may approximate the product VLV in Eq. (1) by the form (/?(r — r') ip{s' — 
s) A(r' — s', r), where '^{r) reflects the symmetry of the short-range correlations 
and A(r, r) describes their strength and long-range properties. Thus Eq. (1) can be 
rewritten in Fourier-Matsubara representation as: 

E(k, = k^TY. f J^ /i^)Mk+ 2- (2) 

(27r)2 -ic^^ + e(q) 



where e(k) is the free dispersion. The self-energy recovers the BCS form when 
A is independent of distance and time, i.e. A oc S{\s. + q)(5t^„ in Eq. (2). In 
this formalism, the superconducting order is therefore related to the long-range 
properties of the function L, rather than to the strength of the superconducting 
correlations. 

Along the lines of Ref. [8], we use the following phenomenological model for A: 

A(r, r) = A2e-^/^«+A?e-'-/«('^), (3) 



where static (r-independent) correlations are assumed for simplicity. The first term 
corresponds to the local correlations of range and the second term describes 
the phase fluctuations regime above — C(^) is the correlation length of the KT 
theory — and the phase coherence below Tc- Both types of correlations are supposed 
to have the same symmetry. With this model, the self-energy on the real-frequency 
axis becomes Il(k, lo) = ip'^{k) J A(k + q)/ + iO+ -|- £:(q)]. In the following, 
we consider Aq and Ai as temperature independent parameters, and we calculate 
numerically the spectral function yl(k, oj) = — :^Im [uj + iO'^ — e(k) — S(k, ui)] ^ 
for a uniform two-dimensional square lattice, a d-wave factor (/^(k) = ^(cosA;^; — 
cos ky), and a dispersion e(k) representative of the Bi2Sr2CaCu208_|_5 (Bi2212) band 
structure [9]. 

The resulting ^(k, co) at the antinodal point of the Fermi surface is shown in Fig. 1 
(upper panels) for several temperatures above Tc. One can see marked differences 
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Fig. 1. Upper panels: Spectral function at the Fermi point near (vr, 0) for several 
temperatures above Tc in the strong (left) and weak (right) pseudogap cases. The 
particle-hole asymmetry in the spectra is due to the underlying electronic structure. 
Lower panels: Calculated ARPES intensity below and above Tc. The model param- 
eters are Tc = 80 K, qq = 10 a, and [A^ + Af]^/^ = 40 meV. The curves are offset 
vertically by equal amounts. The thick curves correspond to T = Tc. 

between the strong (Aq > Ai) and weak (Aq < Ai) pseudogap cases. In the former 
case, the spectral function is dominated by the contribution of the local correlations 
and remains basically temperature independent. In the latter case, the quasiparticle 
peak broadens with increasing temperature, and moves toward the Fermi energy. 
It is worth noting, in this connection, that the gap in the spectral function does 
not relate solely to the strength of the correlations — or to the amplitude of the 
superconducting order parameter — but also to the range of these very correlations; 
this is particularly clear in the right part of Fig. 1, since Ai is the same for all curves. 
Therefore, the spectral function alone does not provide a reliable measure of the 
amplitude of the fluctuating order parameter above Tc. 

In order to compare our results with the experimental data, we must take into ac- 
count the angular (69 ~ 5°) and energy {Suj ~ 15 meV) resolutions of the ARPES 
measurements. We thus estimate the ARPES intensity /(k, cu) as a convolution of 
the calculated A(k, uj) f{uj), f the Fermi function, with a Gaussian of appropriate 
width in k and u>. The calculated /(k, uj) is shown in the lower panels of Fig. 1. 
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Fig. 2. Upper panels: Density of states from Tc to room temperature in the strong 

(left) and weak (right) pseudogap cases. Lower panels: Calculated tunneling con- 
ductance from r = 10 to 300 K. The model parameters are the same as in Fig. 1. 
The curves are offset vertically by equal amounts. The thick curves correspond to 

r = Te. 

Because of the k-dependence of the spectral function, the effect of a finite angu- 
lar resolution is mainly to add a background in the ARPES signal. In the strong 
pseudogap case, wc find that the peak position is more or less temperature indepen- 
dent, in agreement with the experimental results in underdoped Bi2212[2]. In the 
weak pseudogap case, however, the peak approaches the Fermi level as observed in 
overdoped Bi2212[2]. The ARPES lineshapes in Fig. 1 do not show a sudden disap- 
peareance of the quasiparticle peak at Tc [10]. This might be due to our assumption 
of static correlations in Eq. (3). Short time correlations above Tc can be expected 
to decrease the quasiparticle lifetime and eventually to induce a redistribution of 
the spectral weight over a large energy window. 

In Fig. 2 (upper panels), we plot the calculated density of states (DOS) N{io) = 
J d'kA(k, uj) for temperatures above Tc (the DOS below Tc is equal to the DOS at 
Tc in the present model). The lower panels show the tunneling conductance across 
Tc, calculated as dl/dV oc f def'{e — eV) f dk [Tkp ^(k, e) with a uniform matrix 
element = Tq. The behavior of the DOS as a function of temperature parallels the 
behavior of the spectral function in Fig. 1. In the strong pseudogap case, the effect 
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of the local correlations is to broaden the DOS, resulting in smaller coherence peaks 
and in a large zero-bias conductance. In addition, we find characteristic shoulders 
near the Fermi energy, reminiscent of the low-energy states observed in Bi2212 
vortex cores[3,4]. As Fig. 2 shows, these states disappear above Tc, leading to a more 
U-shaped spectrum. Moreover, we have verified that these structures are absent if 
Qo ^ 20 a, and that they do not appear if Aq = or Ai = 0. These findings indicate 
that the low energy states result from the interplay of the short-range and long-range 
correlations. This is confirmed by our numerical results in the weak pseudogap case. 
Indeed, subgap structures can be seen in the DOS for T > 100 K, i.e. a temperature 
range where ^(T) is smaller than ^o- 

In the tunneling spectra shown in Fig. 2 (lower panels), the superconducting gap 
evolves smoothly across Tc into a pseudogap, but the coherence peaks and the gap 
structure is suppressed more rapidly in the weak pseudogap case as the tempera- 
ture is raised. Also, the width of the pseudogap appears to be almost temperature 
independent in both strong and weak pseudogap cases. These trends are consistent 
with the experimental observations in under- and overdoped Bi2212[l], respectively. 

In summary, based on a phenomenological theory of HTS connected to the prop- 
erties of the Cooperon propagator, we have investigated the effect of local pairing 

correlations and long-range phase fluctuations on the spectral properties of a 2- 
dimensional d-wave superconductor. The energy and temperature dependence of 
the spectral functions was shown to vary with the relative strength Aq/Ai and 
the relative range go/^{T) of both types of correlations. In particular, low energy 
states were found to appear in the DOS in two distinct regions of the parameter 
space: Aq/Ai > Qo/S,{T) <^ 1 (strong pseudogap case at low temperature) and 
Aq/Ai < 1) Qo/S,(T) > 1 (weak pseudogap case at high temperature), i.e. when one 
type of correlation is stronger in strength and weaker in range with respect to the 
other type of correlation. 
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